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Abstract. We consider the role of rescattering of secondary kaons on the dilepton branching ratio of the 
4> meson. In-medium mass modifications and broadening of kaons and <j> mesons are taken into account. 
We find in the framework of a Bjorken scenario for the time evolution of the expanding fireball that the <j> 
yield from dimuons is moderately or at least only slightly enhanced compared to that from kaon-antikaon 
pairs. The relation to experimental yields measured by the NA49, NA50 and CERES Collaborations at 
CERN SPS and the PHENIX Collaboration at RHIC is discussed. 

PACS. 25.75.-q Relativistic heavy ion collisions - 14.40.Aq pi and K mesons - 14.40.Ev phi mesons - 
13. 75. Lb meson-meson interactions 



1 Introduction 

Relativistic heavy ion (RHI) collisions present an unique 
opportunity to create nuclear matter at high density and 
temperature where hadron properties become different and 
the phase transition to quark matter with signature of the 
deconfmcmcnt and restoration of chiral symmetry is ex- 
pected. Strangeness enhancement in relativistic nucleus- 
' nucleus collisions compared to nucleon-nucleon collisions 
has been suggested as a possible signal for the formation of 
the quark gluon plasma (QGP) [Tj. The dominant produc- 
tion of ss pairs via gluon-gluon interaction in the plasma 
may result in an enhanced number of strange and multi- 
strange particles produced after hadronization. In partic- 
ular, the production and decay of the <fi meson have been 
long recognised as an important probe for the state of 
matter produced in RHI collisions PEj. The possible en- 
hancement of its yield has been suggested as an important 
signature of the formation of a QGP state 2 , in which free 
ss pairs would coalesce to form <j) mesons, whereas their 
production in pp collisions is suppressed according to the 
Okubo-Zweig-Iizuka rule [3]. 

Phi meson production was measured in central Pb+Pb 
collisions at 158 A GeV at the CERN SPS. The NA49 Col- 
laboration identified the phi meson via the decay channel 
<p -> K+K~ g], while the NA50 Collaboration identified 
it via the 4> ~^ decay 5 . It was found that the ex- 

tracted number of phi mesons from the dimuon channel 
exceeded those number extracted from the K + K~ chan- 
nel by a factor between two and four jS] . 

It has been argued that this difference could be at least 
partially attributed to the fact that not all phi mesons can 
be reconstructed via invariant mass technique applied to 
K + K~ pairs. Due to the rescattering of the K + and K~ 



in the medium, the reconstructed invariant mass could re- 
sult outside the original phi meson peak in such a way 
that some <f> — > K + K~ events contribute to the exper- 
imental background rather than to the measured yields 
[7|. The modification of the visible <j> spectra due to kaon 
rescattering was first investigated within a basic RQMD 
model, showing that kaon rescattering can account for 
the observed rapidity distributions, but not for the dif- 
ferences in the m t spectra, neither the inverse slopes nor 
the relative yields, with a resulting suppression of the <f> 
meson yield via the hadronic KK channel of 40-60% [7|. 
This effect has been further studied in the framework of a 
spherically expanding fireball model by including different 
possible changes (increasing/decreasing) of kaon masses in 
the hadronic medium |H] , as well as in the framework of a 
multi-phase transport (AMPT) model by including den- 
sity dependent in-medium modifications of kaons and <f> 
mesons [2]. In none of these cases the large difference be- 
tween the two SPS experiments could be explained. Such 
a large enhancement of four at low transverse mass from 
dimuon over the dikaon channel could only be explained 
in the AMPT model by an ad hoc increase of the phi 
meson number from the initial stage by about a factor 
of three, that might then suggest that other mechanisms, 
such as formation of colour ropes or quark-gluon plasma, 
are needed for 4> meson production 0. 

To shed light on this so called "</> puzzle", two ex- 
perimental collaborations planned to measure the <fi yield 
simultaneously in both 4> — > e + e~ and KK channels at 
two different energy regimes, namely the CERES Collab- 
oration at CERN SPS and the PHENIX Collaboration 
at the BNL Relativistic Heavy Ion Collider (RHIC). The 
CERES Collaboration measured <j> meson production in 
central Pb + Au collisions at Ei a b = 158 A GeV. Recently 
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reported data show that the yields and inverse slope pa- 
rameters in both decay modes agree within errors. A yield 



the 



decay mode larger than 1.6 times the yield 



in the K+K~ channel is excluded at 95% CL [H]. The 
PHENIX Collaboration measured the <j> yield in Au + Au 
and d + Au collisions at ^/sii = 200 GeV [TTIT2IT5IT31 . 
Whereas definitive results on the <fi yield at mid-rapidity 
in the Au + Au collisions as measured in the K + K~ chan- 
nel have been already reported |12|. only preliminary re- 
sults on the corresponding dilepton channel are available 
so far |11I13I14| . According to recently reported prelim- 
inary production results, the temperatures measured in 
the two decay channels are in agreement within errors, 
while the yield in the e + e~ channel seems to be larger 
than in the K + K~ decay channel. However, the errors in 
the dielectron channel, both statistical and systematic, are 
still too large for a definite statement (Tl]. On the other 
side, recently reported preliminary data on the yields of 
both <j) — > K + K~ and <j> — > e + e~ in d + Au collision show 
that the two data samples are consistent in both the total 
invariant yield and inverse slope |13| . 

In heavy ion reactions at energies reached at RHIC, 
the medium after hadronization is baryon dilute but me- 
son rich |15|. Since by far the most abundant particles 
are pions, one can speak about pion matter created in 
heavy ion collisions at RHIC energies. The kaon-pion gas 
was studied in detail within chiral perturbation theory 
(ChPT) in ref. JSj, focusing thereby on the quark con- 
densates. ChPT, proposed for the description of interac- 
tions of pseudo-scalar mesons at low energies, although 
being an adequate tool for studying the problem at low 
temperatures (T < M„), cannot be applied at typical tem- 
peratures reached at RHIC (T ~ 200 MeV). At high tem- 
peratures, more phenomenological approaches come into 
play and ChPT provides a useful tool to control the low 
temperature limit of such phenomenological models. 

In ref. |17j , the modifications of kaon properties in iso- 
topically symmetric hot pion matter has been calculated. 
Starting from kaon in-medium masses and mean field po- 
tentials calculated in isotopically symmetric pion matter 
to one loop of chiral perturbation theory, the results have 
been extended to RHIC temperatures using experimen- 
tal data on wK scattering phase shifts. Kaon in- medium 
broadening and <f> in-medium width were calculated as 
well. Here we are using the results of ref. ^2] to investi- 
gate the role of the rescattering of secondary kaons inside 
the pion matter on the apparent dilepton branching ratio 
of the 4> meson. 

The paper is organised as follows. In section PJ] we 
briefly review the main results of kaon in-medium masses 
and broadening, mean field potential, as well as of (f> me- 
son in-medium width in an isotopically symmetric hot 
gas of pions investigated in ref. ^Jj. In section we dis- 
cuss the identification of <p mesons from both the dimuon 
and dikaon channels and calculate its apparent dilepton 
branching ratio in the framework of an evolutionary Bjorken 
scenario for the expanding fireball. Summary and conclu- 
sions are finally given in section 01 



2 In-medium modification of kaons in pion 
matter 

We start with the <f> meson in-medium width, which can 
be written in the form 



/ 



*KK 
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contains the in-medium changes of the <p meson decay 
width. Here p* 3 (Af , M K , M K ) and p* 3 (M*,M K , M K ) are 
the kaon three-momenta in the (f> meson rest frame for an 
in-vacuum and in-medium decay respectively. is the 
kaon collision width. The in-medium masses are Ml = 
M^ - 2V K and M* K = M K + SM K , where V K is the mean 
field potential and 8Mk is the mass shift for kaons in iso- 
topically symmetric pion matter |17| . 
The kaon collision width is given by 



r. 



i 



K 



6M^(2tt) 2 ^ 



dn* 



i$ 2 (Vs,M„,M K ) , 



(3) 



where A^(s) are the on-shell partial wave projections of 
the irK amplitudes with total isospin 7=1/2 and 3/2. 
The function $ 2 (yfe,M„,M K ) = Ttp* 3 (-\/s, M n , Mk)/\/s 
is the invariant nK phase space with the cm. momen- 
tum of the irK system p*(y/s, M^, Mk)- The scalar pion 
density is defined by dn S7T — dn V n/(2Ex) via the number 
densities of pions given by the Bose distribution function 



dn V7T — 



<Ppk 

(2tt) 3 



I j, I - 1 



(4) 



-fin- is the 7r + chemical potential and 



where 
/J^o = 0. 

The on-shell ixK amplitudes have been calculated in 
ChPT to the order p 4 by several authors (see e.g. [T3j and 
references therein). To lowest order p 2 , the off-shell wK 
amplitudes are given in ref. |19| . The self energy operator 
£(Mg, Mk), mass-shift 5Mk and mean field potential 
Vk for kaons in an isotopically symmetric pion gas have 
been calculated to leading order in pion density in chi- 
ral perturbation theory from the irK forward scattering 
amplitudes determined at the one loop level. The ampli- 
tudes were expressed in terms of the s- and p-wave ttK 
scattering lengths and the s-wave effective ranges. The 
derived expressions, valid at T < M„, can be extended 
via a more phenomenological approach to typical chemi- 
cal freeze-out temperatures at RHIC (T - 170 MeV PJOj) 
by rewriting the s-wave parts of the amplitudes A in 
terms of phase shifts and by assuming the behaviour of 
the p-wave to be fixed by the a[ scattering length and 
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Fig. 1. Kaon collision width r£ (left) and in-medium <f> meson width r^_ >KR (right, solid line) versus temperature T in 
isotopically symmetric pion matter. The dashed curve in the right panel corresponds to a phenomenologically increased (j> width 
r^ oi = r% + r^+ l p, in which an additional collisional broadening by p mesons is taken into account (see the detailed discussion 
in section |2J . 



the resonance if*. The resulting amplitudes satisfy unitar- 
ity. The experimental 7rif scattering phases are then well 
reproduced; the low-temperature limit matches smoothly 
with one-loop ChPT. For the details, we refer to ref. |17j . 
The resulting kaon collision width and in-medium <p me- 
son width are shown in fig.^ a s a function of temperature 
in the left and right panel respectively. 

3 Apparent dilepton branching ratio of the o 
meson 

Since the <f> meson is unstable, it can only be detected 
from its decay products of either the kaon-antikaon pair 
or the lepton pair. Dileptons from <p decays leave the pion 
matter essentially undistorted by final state interactions, 
whereas the secondary kaons rescatter with the medium 
and can contribute to the experimental background. We 
investigate whether this can result in an increase of the 
apparent dilepton branching ratio or not. 

Let us assume that <j) mesons leave the heavy ion re- 
action at time r. In the simple case of constant widths, 
as already pointed out in ref. |17|. the probability for the 
kaons to leave the reaction zone without rescattering is 
given by 

w ~ / exp (-2 r* K (r - t)) cxp (-JJ t) T dt . (5) 
Jo 

The first factor in this expression is the probability for two 
kaons to escape from the reaction zone without rescatter- 
ing. The second factor is the survival probability of the (j> 
meson at time t. The last factor 1^ dt is the probability to 



decay into the kaon pair during dt. Notice that r£ dt has 
the meaning of a decay probability into the kaons which 
rescatter with pions from the surrounding medium. Se- 
lecting kaon pairs with invariant masses of the <fi meson 
which suffered no rescatterings, the number of (f> mesons 
observed in the K K channel equals 

N KR ~ exp (-r; t) + -ft. [exp r) 

K <fi 

- exp (-2 II r)] . (6) 

The first term arises due to the vacuum decays, whereas 
the second term is given by eq. ©. The number of <j> 
mesons observed in the dilepton channel equals 

N e+e - ~ B exp (-r; t) +B'[1- exp (-2^ r)] , (7) 

where B and B* are the vacuum and in-medium dilepton 
branching ratios of the (f> meson respectively. We assume 
that the dilepton channel suffers no in-medium modifica- 
tions, i.e. B* I B = r^/r^. Notice that the second term of 
eq. 10 can be also written as 

B* [1 - exp (-r; t)] = / exp (-^ t) 2^ e+e - dt 

Jo 




exp (-r; t) b r^dt , (8) 



whose interpretation is straightforward, being exp(— i^ t) 
the survival probability of the <fi meson at time t and 
r$— >e+e~ dt its probability to decay into a dilepton pair 
during dt. 



4 



E. Santini et al: Medium effects on <j> decays to dilepton and kaon-antikaon pairs 



Turning now to temperature dependent, and therefore 
time dependent widths, the previous relations should be 
modified in the following way: 



exp (-2 r£ (r - t)) => exp 
exp (~r£ t) => exp 
exp (-T£t) => exp 



rm dt' 
r;(t') dt 1 
r r;(t>) dt' 



so that the number of <f> mesons decaying into KK pairs 
is determined by 

N KR ~ exp (-£r;(t') dt'^j +jTexp (-2^^) dt'^j 

x exp (-J r;(t') dt'^j r dt , (9) 

and those decaying into dileptons by 



N R+f 



B exp 



rut') dt' 



+ bJ^ exp f- J r;(t') dt'^j r dt . (10) 

Finally, the apparent dilepton branching ratio with respect 
to the vacuum branching ratio is then simply given by the 
ratio of both quantities 10 and l|l(J|) , i. e. 



gapp 

B 



N e+e - 



N 



(11) 



KK 



In central relativistic heavy ion collisions the distri- 
bution of matter is approximatively uniform in rapidity, 
at least in the mid-rapidity region, and the geometry of 
the collision can be considered as cylindrically symmet- 
ric |21| . In the framework of a Bjorken scenario, exten- 
sively used to describe the space-time evolution of ultra- 
relativistic heavy ion collisions, the time and temperature 
of the evolving system are related according to [2] by 



(12) 



where c 2 s is the sound velocity squared for a baryon less 
gas, and To is the initial temperature of the system that is 
assumed to be fully thermalized after the time to- Because 
our investigations are restricted to matter created in heavy 
ion reactions after hadronization, we have applied eq. i|12[) 
to describe the temporal evolution of the hadronized sys- 
tem reinterpreting to and To by the hadronization time 
th and hadronization temperature Th respectively. For an 
ideal relativistic gas, the energy density and the pressure 
are related by e = 3p, therefore one finally has 



- = — = 3 

c 2 s dp 



and 



T(t) 



-1/3 



(13) 



In section |21 expressions for the kaon collision width 
and in-medium (f> meson width as functions of the temper- 
ature were given. The apparent dilepton branching ratio 
can be evaluated from eq. l(TT|l using eqs. © and ijTUI) by 
performing the t —* T variable replacement according to 
the transformation law (|13fl . In addition to the onset of the 
T(t) evolution, the parameters Th and in particular th de- 
termine mainly the slope of the temperature decrease with 
ongoing expansion. Thus, both parameters may strongly 
influence the final results in the following way: the larger 
Th and th are, i.e. the slowlier the hadronic phase cools 
down, the stronger the in-medium effects are. The result- 
ing ratios B app /B for three different sets of these param- 
eters are shown in fig. El as a function of the phi mesons's 
escape time r. First of all, we chose in our model cal- 
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Fig. 2. Ratio between the apparent and vacuum dilepton 
branchings B app /B versus the 4> meson's escape time r (solid 
curves). The dashed curves correspond to the case in which 
an additional collisional broadening in the phi meson width 
is taken into account (see right panel of fig. 0. The verti- 
cally hatched area corresponds to a region around a r value 
of « 32 fm/c predicted by transport calculations for RHIC 
conditions. The horizontally hatched band at B app /B — 1.6 
indicates the upper limit of the enhancement factor according 
to the observations of the CERES Collaboration at SPS (see 
discussion below). 



dilation an initial temperature of Th — 200 MeV and 
a corresponding time of th = 1 fm/c leading to a fast 
cooling behaviour of the expanding fireball. In this quasi 
extreme scenario, we find that the phi meson yield from 
the dilepton channel is only slightly higher by a factor of 
about 1.25 than that from the KK channel - at least in 
the physically interesting r region. According to transport 
calculations, the term exp(— tJ^) is of the order ~ 1/2 at 
RHIC energies *-e. t ~ 32 fm/c. Therefore, the ex- 
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pected region of physical interest is depicted in fig. as 
vertically hatched area. The B app /B ratio from a model 
calculation with practically the same fast cooling, but a 
smaller and therefore more realistic hadronization temper- 
ature of T/j = 170 MeV according to transition temper- 
atures obtained from 2-flavour lattice QCD [23j is even 
essentially smaller, because the interval of the tempera- 
ture integration is shrinked. The consequences of a larger 
hadronization time are also depicted in fig. [21 With a value 
of th — 5 fm/c, the temperature of the expanding hadronic 
fireball decreases much slowlier compared to the previous 
parametrizations. In this case, the in-medium effects influ- 
ence the ratio between the apparent and the vacuum dilep- 
ton branchings B app / B in a much stronger manner, that 
overcompensates the (converse) effect of an even smaller 
hadronization temperature Th — 150 MeV according to 
transition temperatures obtained from 3-flavour lattice 
QCD [22] • The enhancement factor in the latter scenario 
turns out to be about 1.5 — 1.6 in the region around a r 
value of w 32 fm/c. 

Nevertheless, all these values are in qualitative agree- 
ment with the observations of the CERES Collaboration 
at SPS energies, who excludes an enhancement factor larger 
than 1.6 at 95% CL 1 , and contrasts the strong enhance- 
ment factor of 2 — 4 resulting from the comparison between 
the yields measured for the p + p,~ and K + K~ channels by 
the NA50 and NA49 Collaborations respectively. Well, a 
large enhancement of that magnitude is still not totally 
excluded within the simplified view of a Bjorken scenario, 
that assumes only a longitudinally expanding medium. We 
want to emphasise once more that the question, how fast 
the temperature of the hadronic fireball decreases, seems 
to be very crucial. A more realistic description taking ad- 
ditionally a transversal expansion into account might pro- 
mote a fast cooling scenario, that tends to result in not 
so strong enhancement factors. The exclusion of a large 
enhancement due to in-medium effects and kaon colli- 
sional broadening in heavy ion collisions at RHIC ener- 
gies results also from multi-phase transport calculations. 
In particular, the production of </> mesons reconstructed 
from dikaon and dilepton decays has been investigated in 
ref. [Hj for central Au+Au collisions at a RHIC energy of 
•y/sNN = 130 GeV. There have been included density de- 
pendent in-medium modifications of kaon masses derived 
from a chiral effective Lagrangian and a linear decrease of 
the 4> mass according to a scaling law » m^(l— ap/po), 
with a = 0.0255 and the normal nuclear matter density 
Pq. It was found that the <fi meson yield from the dilep- 
ton channel is about 1.5 times larger than that from the 
K + K~ channel in the absence of additional <fi meson pro- 
duction due to the formation of an initial partonic stage. 

The idea of a (f> meson which decouples easily from the 
hadronic fireball is widely diffused, as a consequence of 
the small cross section for its scattering with non-strange 
hadrons. However, available calculations give contradic- 
tory results 24 25 26 ,27 . Whereas some of them support 

1 The cj) yields in the leptonic and the hadronic channel pub- 
lished in ref. HOI result in a ratio B app /B of 1.0 ±0.31 (stat) ± 
0.28 (syst), i.e. a maximum value of about 1.6 (upper limit). 



this idea |24U25U26j . a calculation based on a Hidden Lo- 
cal Symmetry Lagrangian shows that the collision rates of 
the 4> meson in a hot hadronic matter of pseudo-scalar (jr, 
K) and vector mesons (p, w, K* , 0) is rather large [2IJ 
and dominated at temperatures between 150 MeV and 200 
MeV by the K* followed by K and p contributions, while 
the contributions from 7T, to and <f> are smaller. The cor- 
responding <fi mean free path at temperatures above 170 
MeV is between 2.4 fm and 1 fm, and therefore consider- 
ably smaller than the typical size of the hadronic system 
created in heavy ion collisions (~ 10—15 fm). This would 
imply that most <j) mesons produced after the hadroniza- 
tion would not leave the hadronic system without scatter- 
ing. It was found, that reactions changing the number of 
cf> mesons account for more than 80% of the total collision 
rate, so that at the early stage of the expansion, when the 
effect of the decay is negligible, the collision rate is mainly 
responsible for the decrease of the <f> number with respect 
to its value at the hadronization. 

Since a pion gas is only the first approximation to the 
hadronic final state of a relativistic heavy ion collision one 
can estimate the stability of the present investigations by 
the inclusion of the 4> + p — > KK channel, i.e. the next im- 
portant channel which gives - according to the results in 
ref. [23 - rise to in-medium modifications of the (j> — * KK 
decay. However, in contrast to the influence of the pion gas 
|17| we consider the contributions from the p mesons in 
a phenomenological way: following ref. |23 the collisional 
width .T^+p due to the scattering of <j>'s on p mesons is 
added to the in-medium width T£ from the interaction 
with the pions. To do so we extrapolated the results of 
ref. 27 imposing the boundary condition -Tf^J, = at 
T = as well as a fast decrease to zero for T — > and cal- 
culated the apparent dilepton branching ratio using this 
time a total <f> meson width equal to T^ = T5 + r?f . 

Fig. H shows the obtained total cj> meson width r^ ot as a 
function of temperature. Obviously in a very high but nar- 
row temperature range, the + p collision rate dominates 
on the in-medium <f> decay width. When the system cools 
down, the additional rate drops fast and becomes negligi- 
ble compared to r£ in the lower temperature domain. 

Fig. [21 depicts also the corresponding ratios B app /B 
with comparison to their values obtained when the colli- 
sional broadening due to p mesons is neglected. As one 
can easily see, the introduction of a collisional broaden- 
ing does not significantly influence the apparent dilepton 
branching ratio. The increase of the total <f> width reduces 
simultaneously the phi decays into dileptons as well as 
the number of <f> mesons reconstructed via the KK de- 
cay channel. The two reduction factors turn out to be 
comparable in magnitude, so that at the end the ratio 
B app /B is only very slightly additionally enhanced. This 
fact is retained for all of the three parametrizations being 
under consideration. Furthermore, this fact demonstrates 
that for the present purpose the fireball expansion after 
hadronization can safely be described by an expanding 
pion gas. The reason is simply that the integration over 
the entire space-time volume of the expanding system is 



6 



E. Santini et al: Medium effects on <j> decays to dilepton and kaon-antikaon pairs 



dominated by the lower temperatures where the heavier 
mesons are suppressed by their Boltzmann factors. 



4 Summary and conclusions 

We have investigated the role of the rescattering of sec- 
ondary kaons inside pionic matter on the apparent dilep- 
ton branching ratio of the <f> meson in heavy ion collisions 
at RHIC energies, taking into account in-medium mass 
modifications and broadening of kaons and cf> mesons in an 
isotopically symmetric hot pion gas. We have found within 
the framework of a Bjorken scenario for the temporal evo- 
lution of the expanding hadronic fireball that the phi me- 
son yield reconstructed from the dilepton channel is only 
by a factor between about 1.1 and 1.6 higher than that 
reconstructed from the kaon-antikaon channel. We have 
demonstrated that the magnitude of this enhancement de- 
pends also on the question, how slowly the temperature 
of the pionic medium decreases. Nevertheless, for reason- 
able parametrizations of the cooling scenario, our results 
agree quantitatively well with recent observations of the 
CERES Collaboration at CERN SPS energies. The inclu- 
sion of further in-medium effects through the addition of a 
possible collisional broadening in the <f> meson in-medium 
width does not change considerably our results. For heavy 
ion collisions at RHIC energies, only preliminary results 
on the 4> meson yield measured via the two decay channels 
are available up to now. The errors in the dilepton channel 
are unfortunately still too large for a definitive statement. 
Future data from RHIC will help bring more insight into 
the meson propagation in a dense and hot medium. 
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